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Abstract. Non-local boundary conditions for Euclidean quantum gravity are proposed, 
consisting of an integro-differential boundary operator acting on metric perturbations. In 
this case, the operator P on metric perturbations is of Laplace type, subject to non-local 
boundary conditions; by contrast, its adjoint is the sum of a Laplacian and of a singular 
Green operator, subject to local boundary conditions. Self-adjointness of the boundary 
value problem is correctly formulated by looking at Dirichlet-type and Neumann-type 
realizations of the operator P, following recent results in the literature. The set of non- 
local boundary conditions for perturbative modes of the gravitational field is written in 
general form on the Euclidean 4-ball. For a particular choice of the non-local boundary 
operator, explicit formulae for the boundary value problem are obtained in terms of a 
finite number of unknown functions, but subject to some consistency conditions. Among 
the related issues, the problem arises of whether non-local symmetries exist in Euclidean 
quantum gravity. 
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1. Introduction 

The last decade of efforts on the problem of boundary conditions in (one-loop) Euclidean 
quantum gravity has focused on a local formulation, by trying to satisfy the following 
requirements: 

(i) Local nature of the boundary operators [1-5] . 

(ii) Operator on metric perturbations, say P, and ghost operator, say Q, of Laplace type 
[5]. 

(iii) Symmetry, and, possibly, (essential) self-adjointness of the differential operators P 
and Q [4, 5]. 

(iv) Strong ellipticity of the boundary value problems obtained from the operators P and 
Q, with local boundary operators Bi and B2, respectively [5, 6]. 

(v) Gauge- and BRST-invariance of the boundary conditions and/or of the out-in (one- 
loop) amplitude [4-8]. 

At about the same time, in the applications to quantum field theory and quantum grav- 
ity, non-local boundary conditions had been studied mainly for operators of Dirac type 
(see, however, [9]), relying on the early work by Atiyah, Patodi and Singer on spectral 
asymmetry and Riemannian geometry [10]. What is non-local, within that framework, is 
the separation of the spectrum of a first-order elliptic operator (the Dirac operator on the 
boundary) into its positive and negative parts. This leads, in turn, to an unambiguous 
identification of positive- and negative-frequency modes of the (massive or massless) Dirac 
field, and half of them are set to zero on the bounding surface [8, 11-13]. 

On the other hand, non-local boundary conditions for operators of Laplace type had 
already been studied quite intensively in the literature, from at least two points of view: 

(i) The rich mathematical theory of pseudo-differential boundary value problems, where 
both the difi'erential operator P and the boundary operator B may be replaced by integro- 
differential operators [14]. One may consider, for example, the boundary value problem 
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on a bounded open set f2 C 3?"^ with smooth boundary dfl: 

Pu = f infi (1.1) 

Tu = (fi at dfl (1.2) 

where P is the Laplace operator —A (from the point of view of the leading symbol, it is 
more convenient to define the Laplace operator with a negative sign in front of all second 
derivatives) and T is a trace operator, e.g. Tu = ^qu = [ttjan in the Dirichlet case, or 

Tu = 7iw = U;N in the Neumann case, where -^n denotes the inward-pointing normal 

L ' J dO. 

derivative at dVl. The theory we are interested in includes both the system 

A=(^^ :C'^(p)^C°°(n)xC°°{dQ) (1.3) 

and its solution operator, or parametrix, 

=(^R : C°°m X C^idQ) C°°m. (1.4) 

With this notation, R is the Green operator [14] solving the problem (cf (1.2)) 

Pu = f in n and Tu = at dCl (1.5) 
and K is the Poisson operator solving the problem (cf (1.1)) 

= in fi and Tu = ip at dVt. (1.6) 
The Green operator can be expressed in more detail as 

R = Qq + G (1.7) 
where Qq is the pseudo-differential operator (see appendix A) defined by 

Qf = Cn [ r^^^^dx (1.8) 



truncated to fl (this implies extending / by on ^'^/fl, applying Q, and restricting to 
Q), and G is a special term, called a singular Green operator, adapted to the choice of 
boundary conditions. Thus, to get a general calculus, one has to consider systems of the 
form [14] 

for some integers j,k,j',k' such that 

A : c°°{Tiy X c°^{dn)'' c°°{Tiy' x c°°{dn)''' (i.io) 

where P is a pseudo-differential operator on 3?*^, P+ is its truncation to f2, G is a singular 
Green operator acting in fi, T is a trace operator T : Q ^ dfl, is a Poisson operator 
K : dVt — > O, and 5 is a pseudo-differential operator acting on the boundary of VL. 
The trace operators we are interested in can take the form [14] 

To« = 7o« + T> (1.11) 

or, instead [14], 

Tiu = -iiu^SQ-iQU + T[u. (1.12) 

With this notation, one has [14] 



i;Nyu j = 0,l,.... (1.13) 



Moreover, Tq and T[ are integral operators going from Q, to dO,, and the map Sq acts on 
functions on dO.. For example, in population theory, one studies the condition [14] 

w(0) = / u{t)f{t)dt (1.14) 
Jo 

expressing the number w(0) of newborn individuals as a function of the age profile u{t). 
This is a special case of the homogeneous condition, with Q, = 3?+ (cf (1-11)) 

7o« + T^« = (1.15) 
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with [14] 



/•OO 

-T^u= / u(t)f(t)dt. 
Jo 



(1.16) 



(ii) Bose-Einstein condensation models, where integro-difFerential boundary operators lead 
to the existence of bulk states and surface states [15]. More precisely, given the function 

q e Li(3fJ) n L2(3?), one defines [15] 



1 roo 

M^) ^ ^ E / e-^«/«5(y)dy. (1.17) 



The function is, by construction, periodic with period 2ttR, and tends to g as i? tends 
to OO. On considering the region 

BR={x,y:x^ + y^<R''} (1.18) 

one studies the Laplacian acting on square- integr able functions on Br, with non-local 
boundary conditions given by [15] 



U-N 



+ i qnis - s')u{Rcos{s'/R), Rsm{s' /R))ds' = 0. (1.19) 

9Br JqBr 



In polar coordinates, the resulting boundary value problem reads [15] 



+ -— + ^Tr^]=Eu (1.20) 



du 
dr 



Qj.2 J. Qj. J.2 Q^2 

/7T 
qR{R{ip-e))u{R,e)de = o. (1.21) 
-TT 



For example, when the eigenvalue E is positive in equation (1.20), the corresponding 
eigenfunction reads [15] 

Ul,Eir,^) = Jlir^/E)e'^'^ (1.22) 

where Ji is the standard notation for the Bessel function of first kind of order I e Z. 
On denoting by q the Fourier transform of q, and inserting (1.22) into the boundary 



condition (1.21), one finds an equation leading, implicitly, to the knowledge of the positive 
eigenvalues, i.e. 

y/EJ'i{Ry/E) + Ji{Ry/E)q{l/R)\^0. (1.23) 

The solutions which decay rapidly away from the boundary are the surface states, whereas 
the solutions which remain non- negligible are called bulk states [15]. 

In the analysis of pseudo-differential boundary value problems, one studies the heat 
equation for the operator B: 

^ + Bj u{t) = for t > (1.24) 



with initial condition 



u{0) = uo (1.25) 



where B is an operator acting like P+ + G and with a domain defined by the boundary 
condition Tu = 0. Such an operator is called a realization of P. One of the aims of 
functional calculus is to make sense of the exponentiation e~*^ under suitable assumptions 
on B. More precisely, one wants to investigate e~*^ so as to get detailed information both 
on the solutions in terms of their data and on the kernel of the solution operator and its 
trace [14]. For this purpose, a basic tool is the analysis of the resolvent 

Rx = {B-XI)-^ (1.26) 

which makes it possible to study also other functions of B, defined by the Cauchy integral 
formula 

f(B) = ^l^f{X)RxdX (1.27) 

with 7 a curve in the complex plane going around the spectrum of B. 

So far, we have tried to convince the general reader that there are many good reasons 
for studying pseudo-differential boundary value problems with their functional calculus on 
the one hand, and their applications to Euclidean quantum gravity on the other hand. 
Now we can outline the plan of our paper, which is as follows. Section 2, relying on [14], 
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describes how to build the adjoint of a Laplacian, when non-local boundary conditions 
are imposed. This scheme is then applied to the gravitational field in section 3, with a 
particular choice of integro-differential boundary conditions. A mode-by-mode form of 
such boundary conditions is studied in section 4, when the background consists of the 
Euclidean 4-ball. Concluding remarks are presented in section 5, and relevant details are 
given in the appendices. 

2. The adjoint with non-local boundary conditions 

In the case of the gravitational field, inspired by section 1, we consider a scheme where 
the differential operator on metric perturbations remains of Laplace type (as well as the 
ghost operator), whereas the boundary conditions are of integro-differential nature. This 
means that the full boundary operator, say B^^"^, may be expressed as the sum of a local 
operator, say B^^"^, obtained from projectors and first-order differential operators [4, 5], 
and an integral operator, so that the boundary conditions read (see appendix B) 



Bab'^hcdix) 



Bab'^hcdix) 



dM 



+ 



W,^,''\x,x')hcd{x')dV' 



M 



(2.1a) 



dM 



where dV denotes the integration measure over M. We may now decide, following DeWitt 
[16], that unprimed lower-case indices refer to the point x and primed lower-case indices 
refer to the point x' . This leads to 



BAd{x) 



J dM 



Bab'^hcdix) 



idM 



+ 



wJ'^'hc'd'dV 



M 



(2.16) 



dM 



which is the form of the boundary conditions chosen hereafter. 

Since we are concerned, for simplicity, with operators of Laplace type in a flat four- 
dimensional background (all curvature effects result then from the boundary only), it is 
very important for us to understand the effect of integro-differential boundary conditions 
on such a class of operators, motivated by section 1 and bearing in mind Eq. (2.1b). For 
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this purpose, following [14], we remark that, after integration by parts, one finds the Green 
formula for P = —A, u e D{P), and v in the domain D{P*) of the adjoint of P: 



{Pu, v)n = ^ - Aw, vj^ = (u,- A vj^+(llpu, pv 



(2.2) 



where the Green matrix [14] reads, in our case. 



U = i 



. / / 



/ 



(2.3) 



whilst p is the (Cauchy) boundary operator, whose action reduces to 



pu = (7ow,7iw). 



(2.4) 



The same property (2.4) holds for v e D{P*). Suppose now that the boundary conditions 
are expressed in the integro-differential form (cf (1-11)) 



7ow + TqU 



= 0. 



The term (upu,pv^ in Eq. (2.2), which is equal to 

(Upu, pv ) = z(7itt, 7of )afi + ^(7o^^, liv)dQ 



(2.5a) 



can be then re-expressed as 



[Upu.pv] = z(7iw,7of)an + *(-ToW, 7if)an 

\ / do, 



(2.56) 



which implies that P*, the (formal) adjoint of P, can be obtained by adding to —A a 
singular Green operator, i.e. [14] 



(2.6) 



supplemented by the local boundary condition 



7ov = at dO,. 



(2.7) 
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By contrast, if the boundary conditions (cf (1-12)) 



7iM + So^ou + T[u 



= 

90 



are imposed, which modify the standard Neumann case, it is convenient to re-express 7itt, 
at the boundary, in the form 

7itt = -So-iQU - T[u (2.8) 
and insert Eq. (2.8) into Eq. (2.5a). This imphes that the adjoint of P now reads 

P*v = - Av + iTi*-fov (2.9) 

subject to the local boundary condition 

7iv = at dn. (2.10) 

In other words, we are discovering a property which is known to some mathematicians, 
but not so familiar to physicists: if an elliptic differential operator (here taken to be 
of Laplace type) is studied with integro-differential boundary conditions, its adjoint is a 
pseudo-differential operator, subject to local boundary conditions. 

Self-adjointness problems are properly formulated by studying the realization of the 
operator P [14]. In our case, this means adding to the Laplacian a singular Green oper- 
ator, and considering a trace operator which expresses the integro-differential boundary 
conditions. More precisely, a Dirichlet-type realization of P = — A is the operator [14] 

Bd=(- A + Gd) (2.11) 

where (see (1-11)) 

Gd = i^i7i + G' (2.12) 

To = 7o + T^. (2.13) 

In our paper, the Poisson operators Ki, for i = 0, 1, are completely determined by the 
requirement of self-adjointness. Indeed, the domains of Bd and its adjoint coincide [14] if 
and only if (cf (2.6) and (2.12)) 

Ki = i T^* (2.14) 
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G' = G'* 



(2.15) 



Moreover, a Neumann-type realization of P = —A is the operator 



B 



n=[-A + Gn) 



Ti 



where (see (1-12)) 



Gn = Ko^o + F' 
Ti = 71 + ^070 + 



(2.16) 

(2.17) 
(2.18) 



FoUowing [14], we use the notation (2.11) and (2.16) for particular realizations of the 
Laplacian, but a notation along the lines of (1.9) if we want to stress the properties of a 
system in the functional calculus of the boundary value problem. The domains of and 
its adjoint are found to coincide [14] if and only if (cf (2.9) and (2.17)) 



Ko = i Ti* 

So = —Sq 
pi _ pi* 



(2.19) 

(2.20) 
(2.21) 



3. Application to the gravitational field 

In the case of the gravitational field, our boundary operator (2.1b) corresponds to the 
trace operator (2.13). The local boundary operator B^^^'^ is taken to be the one for which 
the following conditions are imposed on metric perturbations on a 3-sphere boundary of 
radius a [2]: 



h 



= 



J dM 



Oz 



dM 



= 



dhoQ 6 9 ^ 

or T or 



= 



(3.1) 



(3.2) 



(3.3) 



dM 
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where r e [0, a]. Equations (3.1)-(3.3) express, to our knowledge, the only set of local 
boundary conditions which are of Dirichlet type on hij and hoi, and for which strong 
ellipticity of the boundary value problem is not violated [5, 6]. Since we only want to 
modify the Dirichlet sector of such boundary conditions, which is expressed by (3.1) and 
(3.2), we have to require that (see (2.1b)) 



W, 



c'd' 



00 



Vc',d' 



(3.4) 



Thus, we eventually consider the system (cf our Eq. (1.9), and Eq. (1.6.84) in [14]) 



A 



-A+G\ 
T J 



(3.5) 



where T is of the type (2.13) in its ij and Oi components, i.e. 



hij{x) 



J dM 



+ 



M 



dM 



(3.6) 



dM 



h 



Oi 



dM 



+ 



Wo/'''h,,d'dV' 



M 



dM 



(3.7) 



and of the type (3.3) (cf Eq. (2.18) and set T{ = therein) in its normal component hoo, 
i.e. 



Too'^hcdix) 



AdM 



dhpo 6 
dr T 



u d ■■ 



(3.8) 



dM 



Moreover, —A is the standard Laplacian on metric perturbations in flat Euclidean 4-space, 
and G may be viewed as the direct sum of and Gn (cf example 1.6.16 in [14]), with 



Gd = K^li + G' 



Gn = F' 



(3.9) 
(3.10) 



subject to the self-adjointness conditions 



{K^)Ad 



— I 



J dM 



M 



W^/'^'hc'd'dV 



(3.11) 



dM 
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(i^l)o/'/^cd 



[ Wo/'''h,,d'dV' 
Jm 



(3.12) 



dM 



dM 



{Ki)oo'"^ hcd 
n' ^'^ — ( n' ^'^\ 

^ ah — \y ab J 



= 



(3.13) 
(3.14) 

F'ab'={F'ab'y- (3.15) 



Note that, by virtue of Eq. (3.4), the counterpart of T{ (see (2.18)) vanishes in our problem 
(as we said after (3.7)), and hence Kq vanishes as well (see (2.19)), so that Gn reduces 
to F' as we write in (3.10). Furthermore, the condition (2.20) is satisfied by virtue of the 
boundary condition (3.8). 

4. Mode-by-mode equations on the Euclidean 4-ball 

At this stage it can be helpful to write down a set of equations for perturbative modes 
of the gravitational field, once that the right-hand sides of (3.6) and (3.7), and the left- 
hand side of (3.3), are set to zero at the boundary. As in section 3, our background is a 
portion of flat Euclidean 4-space bounded by a 3-sphere (hereafter, x corresponds to local 
coordinates on the 3-sphere; x and r, altogether, correspond to the symbol x used before). 
This is quite important in (one- loop) quantum cosmology [7, 8, 11-13], and also as a first 
step towards more complicated field-theoretical models. Under the above assumptions, the 
expansion of metric perturbations on a family of 3-spheres centred on the origin reads [8] 



/ioo(^,T) = 5]a„(T)Q(-)(f) 



n=l 



hoi{x,T) = 



n=2 



bn{r)^^^+cr,ir)si-\x) 



(4.1) 



(4.2) 
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hij{x,T) = ^Unir) 
n=3 

oo ^ 



n=l 
oo 



(4.3) 



n=3 



Here, with a standard notation, Q^'^\x), sl^\x) and G\^\x) are the scalar, transverse 



vector, transverse-traceless tensor hyperspherical harmonics on a unit 3-sphere (with metric 
Cij), respectively. On denoting by a parameter with dimension [length] by Ir the 
modified Bessel function of first kind and order r, and on choosing the de Bonder gauge- 
averaging functional, one then finds for the transverse-traceless modes the formula [8] 



Znir) = anTlnilJ^T) 



whilst the vector modes read [8] 



C2(t) = eh{tJiT) 

,n-^n— 1 



fnir) = T 



ei,r 



(n + 2) 

and the scalar modes are given by [8] 

1 r 



--2,n 



(n-2) 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



ei(T) = T SAihif^r) - ^4/3 (iwr)] 
a2(T) = ^ Bihiiir) + BMl^^)] 



b2{T) = B2l2{l^T) - B4h{l^T) 



62 (r) = T WihifiT) - 2B2l2{fiT) 



(4.8) 

(4.9) 

(4.10) 
(4.11) 
(4.12) 
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anir) = - 2 {PT ) + P4:,nIn+2 {P'T ) 

K{r) = P2,nln{pr) + {n+ 1) p3^nIn-2{P'r) - {u - l)p4,n/n+2 (jWt) 



Unir) = T 

(n + 2) 



-P2,nln{pr) + — ^P3^riln-2{pr) 



(n-2)' 



P4,nIn+2 



(4.13) 
(4.14) 



(4.15) 



en{r) = r ^pl,nIn{P'T) - 2p2,nln{pr) - p3,nln-2{pr) 



- PA,nIn+2{pr) 



(4.16) 



Now we follow the procedure outlined in the introduction in a simpler case (see (1.20)- 
(1.23)), i.e. we insert the mode solutions (4.4)-(4.16) of the eigenvalue equation for metric 
perturbations into the mode- by- mode form of the boundary condition resulting from (3.6), 
(3.7) and (3.3). For this purpose, it is convenient to define 



K0i= / W^oi"''''/ic'd' dV 



M 



Kij — I 

JM 



(4.17) 



(4.18) 



The boundary conditions are then expressed by (3.3) jointly with the equations 



= 



hij ~\~ Kij 



0. 



(4.19) 



(4.20) 



The tensor fields kqi and Kij, representing the non-local contribution to the boundary 
conditions for a symmetric rank-two tensor field, are themselves symmetric. They can be 
therefore expanded on a family of 3-spheres centred on the origin according to formulae en- 
tirely analogous to (4.2) and (4.3), with the modes Cn} replaced by the modes |6n, Cn|, 
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say, and the modes {un, e^, fm ^n} replaced by the modes l^t^, e„, /„, The boundary 

conditions (3.3), (4.19) and (4.20) lead, therefore, to the following set of equations for 
perturbative modes, for all n > 3: 



Ik 



+ -ar 



1 den 



(r = a) = 



(4.21) 

(4.22) 
(4.23) 
(4.24) 
(4.25) 
(4.26) 
(4.27) 

The finite-dimensional spaces corresponding to the modes (4.5) and (4.8)~(4.12) should, 
of course, be treated separately, by requiring that 



[K + K] {t = 


a) 


= 


[Cn + Cn] (r = 


a) 


= 


[Un + Un] {r 


a) 


= 


' ..... + en] (r = 


a) 


= 


[fn + fn]{r = 


a) 


= 


[Zn + Zn] (r = 


a) 


= 0. 



dak 6 1 dck 
dr T ar 



(r = a) = if A; = 1,2 



[b2 + b2]{r = a) = [c2 + d2]{r = a) = 
[e2 + e2](T = a) = 0. 



(4.28) 

(4.29) 
(4.30) 



It is now convenient to write more explicitly the tensor fields Koi and k^ , by using the 
definitions (4.17) and (4.18) on the one hand, and the expansions analogous to (4.2) and 
(4.3) on the other hand. This leads to 

Koiix,T)= [ IWo/'^'ho'o' + 2WQ,^°''''ho'k' + Wo,'''^'hk'i']dV' 



= E 

n=2 



hn{r) 



[v? - 1) 



(4.31) 
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M 



dV 



n=3 



(n2 - 1) 



n=l 



oo 



(4.32) 



n=3 



One should now insert the expansions (4.1)-(4.3) into the integrals occurring in (4.31) 
and (4.32), hence reading out the modes bn,Cn,Un,en, fn,Zn for a given form of W^^,*^ , 
chosen to be compatible with (3.4), (3.11) and (3.12). Last, such a solution should be 
inserted into the boundary conditions (4.21)-(4.30), which would be then re-expressed, 
in non-local form, uniquely in terms of modes for metric perturbations, which are known 
from (4.4)-(4.16). 

For example, if one assumes that W^j,*^ has distributional nature, and that suitable 
functions {/i, ...,/6} exist (they should be such that the integrals we are going to build 
exist as Lebesgue or Riemann integrals) for which (it is convenient to factorize t'~^ in our 
ansatz, because it cancels exactly a term r'^ from the integration measure dV) 



(4.33) 



(O'fc') _ /-3 



(4.34) 
(4.35) 



where Vi coincides with the covariant derivative on the boundary denoted by |i so far, the 
method described after Eq. (4.32) leads to, for all n > 2, 



>>n(T) 



^ = /°[A(r,r')a„(r') + /.(x,r')AM 



(n2 

+ MT,T')en{T') 



dr' 



(4.36) 
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Cn{r)= [ f2{T,T')Cn{T')dT' 

Jo 



whereas, if n = 1, one finds 



h{T,T')a^{T') + fsiT,T')e,{T') = 



which imphes the non-trivial property 



/i(r,rO ei(T') 



/3(t,t') ai(TO 



Moreover, if one assumes that 



/\r k 



■f\X k s I 
j) 



the same method imphes, for the first few modes. 



e,{T) = £[3U{T,T')a,{T') + fe{T,T')e^{T')\dT' 



MT,T')b2{T')dT' = 



r Mr,T')c2{T')dT' = 

Jo 

and hence, up to a zero-measure set (see (4.5) and (4.11)), 



/5(t,t') = 



and, for aU n > 3, 



e2(T) = £ [3/4(T,T')a2(T') + fe{T,T')e2{T')\dT' 
Un{r) ^ / f6{T,T')Un{T')dT' 

Jo 
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/ 

^0 



a ^ 



en{r) 



3/4(t, T')an{r') + feir, T')en{T') dr' 



(4.49) 



I 



a 



fn{r) 



(4.50) 



Zn{r) 



r h{T,T')Zn{T')dT' 



(4.51) 



The formulae (4.36), (4.37) and (4.48)-(4.51) for the coupled modes should be inserted into 
the boundary conditions (4.22)-(4.27), hence leading to the non-local form of the boundary 
conditions resulting from the assumptions (4.33)-(4.35) and (4.40)-(4.42). Nothing more 
explicit can be said unless one determines the form of all functions {/i, /e}. 

5. Concluding remarks 

Motivated by the recent developments in Euclidean quantum gravity on manifolds with 
boundary [1-9] and Bosc-Einstcin condensation models on the one hand [15], and by the 
progress in the functional calculus of pseudo-differential boundary value problems on the 
other hand [14], we have considered the mixed boundary conditions (3.3), (4.19) and (4.20) 
for the quantized gravitational field on the Euclidean 4-ball. The main drawback of such 
boundary conditions is, possibly, the apparent lack of an invariance principle leading to 
their derivation. On the other hand, it has been recently proved in [5, 6] that precisely the 
boundary conditions which are completely invariant under infinitesimal diffeomorphisms 
on metric perturbations lead to serious technical problems. In other words, on choosing the 
operator P on metric perturbations to be of Laplace type, the resulting boundary value 
problem fails to be strongly elliptic [5, 6]. This is a technical condition which requires 
that a unique solution should exist of the eigenvalue equation for the leading symbol of P, 
subject to a decay condition at infinite geodesic distance from the boundary and to the 
boundary conditions of the problem [5, 6]. If there is lack of strong ellipticity, the fibre 
trace of the heat- kernel diagonal acquires a non-integrable part near the boundary [5] , and 
hence not even the one-loop semiclassical approximation is well defined (except, possibly, 
if one works with a "smeared" form [13] of the heat kernel for the operator P). 
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In the light of the above difficulties, the possible lack of an invariance principle does not 
seem a sufficient reason for not considering the many interesting motivations leading to the 
boundary conditions (3.3), (4.19) and (4.20). In particular, the Eqs. (4.4)-(4.16), jointly 
with (4.21)-(4.32), supplemented by the self-adjointness conditions (3.11) and (3.12), lead 
to what seems to be a very interesting calculational scheme for quantum gravity on the 
Euclidean 4-ball (see (4.33)-(4.51)). At least six outstanding problems are now in sight: 

(i) Is there a set of non-local symmetries in Euclidean quantum gravity leading to our 
non-local set of boundary conditions (3.3), (4.19) and (4.20) or to a suitable modification 
of our scheme? 

(ii) Is the resulting class of boundary value problems compatible with strong ellipticity? 

(iii) Can one build explicitly a class of bulk and surface states in Euclidean quantum 
gravity with non-local boundary conditions, inspired by the work in [15]? The idea is then 
to prove that, upon inserting (4.4)-(4.16) into (4.21)-(4.32), solutions formally analogous 
to the bulk and surface states of [15] do actually exist (cf (1.20)-(1.23)). 

(iv) Can one make assumptions for W^j^'^ '^' in the boundary conditions (2.1b) which admit 
(4.33)-(4.35) and (4.40)-(4.42) as a particular case? This would elucidate the general 
structure of the kernel of the trace operator in the system (3.5). 

(v) Can one study heat-kernel asymptotics with non-local boundary conditions for the 
gravitational field? 

(vi) Can one perform a path- integral quantization, if the non-local choice (2.1b) is made 
for the boundary data? What form of boundary conditions should be imposed on ghost 
fields? 

If one were able to solve all such problems, an entirely new vision would emerge 
in Euclidean quantum gravity, with a non-trivial impact also on the other branches of 
Euclidean field theories. Thus, encouraging evidence exists that the Euclidean approach 
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continues to play a key role on the way towards further progress in the theory of the 
quantized gravitational field [8, 17]. 
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Appendix A 

To be self-contained, it is important to describe the main ideas behind the definition of 
pseudo-differential operators. For this purpose, following [14], let us recall that the action 
of a differential operator of order m on 3?", 

A{x,D,)= J2 (^aix)D^ (AA) 

\a\<m 

can be expressed, with the help of Fourier transform, as 

where a : (a;,^) — > a{x,^) is the function called symbol or characteristic polynomial, whose 
action is defined by 

a{x,0^ Yl ««(x)r- (^-3) 

\a\<m 

Pseudo-differential operators are obtained by considering, instead of a symbol of polyno- 
mial nature as in (A. 3), a more general function. More precisely, a pseudo-differential 
operator P with symbol p{x, ^) is the operator defined by [14] 

iPu){x) ^ (27r)- / e<--y>h{x, i)u{y)dyd^. {AA) 

Strictly, the definition (A. 4) holds for u e S{^^), but can be extended to a more gen- 
eral class of functions, provided that the symbol p{x^^) satisfies suitable conditions. In 
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particular, it is important to consider symbols which, for some 5 e 3?, are C°° functions 
satisfying the inequality 

\D^^D^p{x,0\<C^A^){l + \^fy-^ ya,p (A5) 

where Cq, /3 is a continuous function. In several applications, one needs also the asymptotic 
expansion of the symbol. For this purpose, one can assume that p{x, ^) is polyhomogeneous, 
i.e. it satisfies (A. 5) and has an asymptotic expansion [14] 

where each p§-i is a C°° function homogeneous of degree 5 — Z in ^ for |^| > 1, and 

P-^P5-i 

is a C°° function satisfying the inequality (A. 5) with 5 replaced by 5 — j, for all j e jV. 
Appendix B 

The construction of the local part of the boundary operator in Eq. (2.1b) deserves further 
comments. In other words, a rigorous formulation of local boundary conditions for opera- 
tors of Laplace type involves vector bundles over M and its boundary with their sections, 
a matrix consisting of projectors and first-order differential operators, and the boundary 
data. A concise description is as follows. We consider an m-dimensional Riemannian man- 
ifold, say (M, ^f), a vector bundle V over M, with a connection V, and operators of Laplace 
type, i.e. 

P = -sr^^'VaVb - E (B.l) 

with E an endomorphism of V. The operator P maps smooth sections of V, say <f, into 
smooth sections of V. In the case of local boundary conditions for P, their general form 
is [5] 



n \ f [cp]dM 

A I-UJ \[if,N]dM 



= {B.2) 
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where 11 is a self-adjoint projection operator, and A is a tangential differential operator on 
the boundary of M: 



A = (i-n) 



(i-n). (S.3) 



With the notation used in (B.3), V is the induced connection on dM, are endomorphism- 
valued vector fields on the boundary, and E is an endomorphism of the vector bundle over 
dM which is a copy of [F]aM, with sections given by [ip\dM- and E are anti-self-adjoint 
and self-adjoint, respectively, and are annihilated by 11 on the left and on the right, i.e. 

n P = P n = (S.4) 

nE = En = o. (S.5) 

As is shown in [1, 4, 5, 7, 8], one arrives at such boundary conditions whenever one tries 
to obtain gauge- and BRST-invariant boundary conditions in quantum field theory. 

In our paper, however, to avoid losing strong ellipticity of the boundary value problem 
in Euclidean quantum gravity [5, 6], we always assume that 

r = (S.6) 

for all « = 1, m — 1. 
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